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SUMMARY 

A previous analysis for fully developed turbulent heat transfer and 
flow with variable fluid properties is extended and applied to the 
entrance regions of smooth tubes and parallel flat plates. Integral 
heat-transfer and momentum eq.uations are used for calculating the thick- 
nesses of the thermal and flow boundary layers. The effect of variable 
properties is deteimined for the case of uniform heat flux^ uniform 
initial temperature distribution^ and fully developed velocity distri- 
bution. A number of other cases in which the fluid properties are con- 
sidered constant are analyzed. The predicted Nusselt numbers for air 
with a uniform wall temperature and uniform Initial temperature and 
velocity distributions agree closely with experimentally determined 
values. 


INTRODUCTION 

Flow and heat transfer in the entrance regions of channels have 
been subjects of considerable interest. The laminar flow case has 
been analyzed in references 1 (pp. 299-310), 2, and 3 (pp. 451-464), 
for instance, and the turbulent case, in references 4 and 5. In these 
analyses constant fluid properties were assumed, and uniform wall 
temperature was postulated in the cases in which heat transfer was con- 
sidered. In the turbulent flow cases the flow was assumed to be tur- 
bulent at all points along the passage. 

The analysis in reference 4 is for a Prandtl number of 1 and is 
based on an assumed l/7-power velocity profile eind the Blasius resis- 
tance formula. The relative changes in heat-transfer coefficient along 
the tube predicted in that analysis are approximately correct, but the 
absolute values of the heat -transfer coefficient at the end of the 
entrance region are about 20 percent too low. The analysis given in 
reference 5 is applicable to cases in which the heat transferred by 
turbulence can be neglected, as might occur at very low Prandtl numbers. 
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The analysis given herein, vhich was made at the NACA Lewis labora- 
tory, is an extension of the analysis given in reference 6 for fully 
developed heat transfer and flow with variable fluid properties to 
the entrance regions of smooth passages. In calculating the growth of- 
the thermal and flow boundary layers in the entrance region, Integral 
eq.uations for heat transfer and momentum are used. For obtaining the 
velocity and temperature distributions in the boundary layers, the same 
assumptions are made for solving the turbulent transfer equations 
as were made in reference 6 . 

The flow is assumed to be turbulent at all points along the 
passage, as in previous analyses, so that the analysis should be 
applicable when disturbances occurring at the entrance are sufficient 
to produce a turbulent boundary layer. According to Prandtl's assump- 
tion, which is in reasonable agreement with experiment, in the presence 
of a laminar boundary layer near the entrance, the turbulent portion of 
the boundary layer behaves as though the boundary layer were turbulent 
all the way from the entrance (ref. 7, p. 74). The present calculations 
may therefore be applicable to the turbulent portion of the boundary 
layer even when a laminar boundary layer exists near the entrance . 


The cases considered in the present analysis are given in the 
following table: 


Phenomenon 

Wall 

Initial 

Initial 

Passage 

Prandtl 

Proper- 

boundary 

temperature 

velocity 


number 

ties 


condition 

distribution 

distribution 




Heat 

Uniform 

Uniform 

Fully 

Tube 

0.73 

Variable 

transfer 

heat flux 


developed 



Variable 

Heat 

Uniform 

Uniform 

Fully 

Parallel 

.73 

transfer 

heat flux 


developed 

flat 

plates 



Heat 

Uniform 

Uniform 

Fully 

Tube 

.73 

Constant 

transfer 

wall tern- 


developed 




Heat 

perature 

Uniform 

Uniform 

Uniform 

Tube 

.73 

Constant 

transfer 

heat flux 





Constant 

Heat 

Uniform 

Uniform 

Uniform 

Parallel 

.73 

transfer 

heat flux: 



flat 

plates 



Heat 

Uniform 

Uniform 

Uniform 

Tube 

.73 

Constant 

transfer 

wall tem- 






Heat 

perature 

Uniform 

Uniform 

Fully 

Tube 

.01 

Constant 

transfer 

Friction 

heat flux 


developed 

Uniform 

Uniform 

Tube 


Constant 



Parallel 

— 

Constant 

Friction 



flat 

plates 
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Although most of the cases in the preceding table are for constant fluid 
properties; the greater portion of the report is concerned with the 
first two caseS; in which the properties are variable* 


MALYSIS 

For calculating heat transfer and friction in the entrance region 
of ducts ; the usual boundary layer assumptions are made; that iS; it is 
assumed that the effects of heat transfer and friction are confined to 
fluid layers close to the surface (thermal and flow boundary layers, 
respectively). The temperature and velocity distributions outside the 
boundary layers are assumed uniform, and the total temperature and total 
pressure are constant along the length of the duct for the region outside 
the boundary layers. More exact analyses (ref. 2) indicate that these 
assumptions are valid, even for laminar flow, except in the region at a 
distance from the entrance .where the boundary layer fills a large portion 
of the tube. For that region, however, the Nusselt numbers and friction 
factors have values very close to the values for fully developed flow. 

In the following analysis, flow in round tubes and between parallel 
flat plates will be considered. The effect of frictional heating will 
be neglected and then later investigated to some extent in the section 
Effect of frictional heating on entrance length. 


Velocity and Temperature Distributions in Boundary i^yers for. Gases 

For obtaining the velocities and temperatures in the flow and 
thermal boundary layers, the differential eq^uations for shear stress 
and heat transfer can be written as follows : 


T 



pe 


du 


d = - k - PCpge^ 

(The symbols used in this report are defined in appendix A.) In the 
preceding equations e and are the eddy diffuslvities for momentum 

and heat transfer, respectively, the values for which are dependent on 
the amount and kind of turbulent mixing at a point. When written in 
dimensionless form, these equations become 

ju du^ P £ du'^ 

^0 ” ^"0 ^ ^0 dy "^ 
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and 


A. - fl- + iL ££_ \ ^ 

<10 " V^o ^0 ^ Pq °Pq ^o/Poj dy"^ 

vhere the Buhecript 0 refers to values at the wall. 

Assumptions . - The following assumptions are made In the use of the 

previous equations for obtaining velocity and temperature distributions 
In the flow and thermal boundary layers: 

1. The eddy dlffuslvitles for momentum e and heat transfer ej^ 
are equal. Previous analyses for fully developed flow In tubes based on 
this assumption yielded heat- transfer coefficients and friction factors 
that agree with experiment for Reynolds numbers above 15^000 (ref. 6), 

At low Reynolds or Peclet numbers (Pe = ReXPr), the ratio of eddy 
dlffuslvitles appears to be a function of Peclet number (ref. 8); but 
for Reynolds numbers above 15,000, It is nearly constant for gases, at 
least for fully developed boundary layers. The same assumption is made 
here for developing boundary layers. 

2. The eddy diffusivity t is given by 

c = n^uy 

in the region close to the wall (y'''<26), and by the Karm& relation 

_ ^2 (du/dy)^_ 

(d^u/dy^) 

in the region at a distance from the wall (y"''>26). These expressions 
have been experimentally verified for fully developed boundary layers 
with variable properties in reference 6. It is assumed here that they 
apply also to developing boundary layers. The quantities n and x are 
experimental constants having the values 0.109 and 0.36, respectively. 

3. The variations across the flow and thermal boundary layers of 
the shear stress x and the heat transfer per unit area q have a 
negligible effect on the velocity and temperature distributions. It is 
shown in figure 12 of reference 6 that the assumption of a linear varia- 
tion of shear stress and heat transfer across the boundary layers 

(x or q = 0 at the edge of the thermal or flow boundary layer) gives 
very nearly the same velocity and temperature profiles as those obtained 
by assuming uniform shear stress and heat transfer across the boundary 
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layers for values of or "between 500 and 5000.^ For small 

values of } such as occur very near the entrance, the effect of 

variable heat transfer will he checked in appendix B and in the section 
Effect of variation of heat transfer across boundary layer on Nusselt 
number for small and constant fluid properties. 

4. The molecular shear stress and heat-transfer terms in the 
eq^uations can be neglected in the region at a distance from the wall 
(ref. 6, fig. 14). 

5. The static pressure can be considered constant across the 
passage. 

6. The Prandtl number and specific heat can be considered constant 
with temperature variation. The variations with temperature of the 
specific heat and Prandtl number of gases are of a lower order of magni- 
tude than the variations of viscosity, thermal conductivity, and 
density. The variations of viscosity and thermal conductivity with 
temperature are assumed to be given by the equations 

^ k/kg = (t/t^)^*^^ 


From the definitions of the quantities Involved it can be shown that 

t/to = 1 - 


Solution of equations . - The details of the solution of the equa- 
tions for shear stress and heat transfer \mder the foregoing assumptions 
are given in reference 6. The resulting relations between t*^ and 


and u“^ and are given in figures 1 and 2. Positive values of 


J 

P 


correspond to heat addition to the gas; negative values, to heat 
extraction. For the cases in which the fluid properties are considered 
constant, the curves for 3 = 0 are used.^ The plotted values of t'*' 
and u"^ eire to be used in the thermal and flow boundary layers, that 


is, for values of j 


■< V 


or y'*‘< 8+. For values of y'*' greater 


^Figure 12 of reference 6 applies to developing "boundary layers If 
the symhols rQ"^ In the figure are replaced "by S'*" or 

^The case In which 3=0 Is a limiting one which can "be approached 
as closely as desired "by making the temperature difference small. In 
the limit, where the temperature difference Is zero or Infinitesimal, the 
properties must be constant. 
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than these values, t^ and are constant and have the values at 

and The relations between and X/D and 6'*' and X/D will 

be obtained in the following sections. 

Developmont of* Thenaal Boundary Layors 

Heat -flow equations for thermal boundary layers . - The heat-flow 
equation for the thermal boundary layer in a round tube can be obtained 
from the following diagram: 



Energy flows into the differential annulus by convection through plane 1 
and leaves through plane 2. In addition to this energy there is a radial 
flow of energy (heat) from the tube surface by conduction and a radial 
flow by convection at radial heat flow due to temperature 

gradient takes place at because the temperature gradient is, by 

definition of the thermal boundary layer, zero at the edge of the 
boundary layer. 

Equating the heat energy entering the annulus to that leaving 
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where the subscripts on the parentheses refer to planes 1 and 2 in the 
diagram; or 



where the differentials of the integrals indicate changes in the X-direction. 
Equation (l) is the desired equation for flow in tubes relating the thermal- 
boundary-layer thickness to the distance along the tube. 


An equation corresponding to equation (l) can be easily derived for 
flow between flat plates. This equation is 


qg dX = d 




( 2 ) 


Equation (2) in the incompressible form has been used extensively 
for calculating the growth of the thermal boundary layer on a flat 
plate (ref, 7; p. 86). 


Uniform heat flux, variable fluid properties (gases) . - The total 
temperature outside the thermal boundary layer does not vary with X 
Inasmuch as no heat penetrates the region outside the boundary layer. 
(The differences between total and static temperature are small except 
at high Mach numbers.) For uniform heat flux at the wall, equation (l) 
can then be integrated to give 


<lQ ^0 X 
°p g 



t§) Pu(rQ - y)dy 


Equation (3) can be written in dimensionless form as 





(3) 


0 


- y'^)dy'^ 


( 4 ) 
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where, from the definitions of the various quantities involved and the 
perfect gas law (constant pressure), 

p/pQ = 1/(1 - (5) 


When equation (4) is to he used , values of and P must first he 

fixed. Values of x/D can then he calculated for various values of 

The relations between t"*^ and and u"*^ and are obtained 

from figures 1 and 2 . 


The equation for flat plates corresponding to equation (4) is 
obtained from equation (2) and is found to be 

p5h^ 


X 

D 



vJo 


(tg^ - t+) 



dy’^ 


(6) 


where D is the equivalent diameter for flow between flat plates, or 
twice the plate spacing. 

Uniform wall temperature, constant fluid properties . - For flow in 
tubes, equation (l) can be written as 



(t - tg) Pu(ro - y)dy 


(7) 


where the quantity in brackets is also the upper limit on the integral. 
For constant fluid properties and uniform wall temperature, equation (7) 
can be written in dimensionless form as 


rc] 



— 1 ^ 

Jo 

^0^ Jo 


where the shear stress Tq Is allowed to vary with X. But 

” ^iAq = 1 - P tg 


since tg does not vary with X. Therefore, 

P = (l/tg+)(l-ti/to) 
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Substituting this expression for P into equation (8) results in 



Equation (9) gives the relation between X/D and 5^"^ for round tubes 

with uniform wall temperature and constant fluid properties. Its use 
is similar to that of equation (4) . 


ITusselt numbers and Reynolds numbers . - It can easily be shown 
(see ref. 9) from the definitions of the quantities involved that the 
Nusselt number and Reynolds number (properties variable or constant) 
for circular tubes with the fluid properties evaluated at the wall tem- 
perature are given by 


NUq = 



( 10 ) 


and 


Rsq = 2u^"^ Tq^ 


( 11 ) 


where 


and 



unro+-y+)dy^ 
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which, by combination with equation (12), becomes 


Ui* = (1-PV) 


05 


^ u^(ro^-y^) ^ ^ 1 

; — dy + 


1-pt' 


l-Ptg 


Orr 


u''’(rQ+-y'^)dy+l 


(13) 


The foregoing equations are for a circular tube. The corresponding 
equations for flow between flat plates are 


RUq 



(14) 


and 



When values of Nuq and I^Sq are calculated, values of rQ and P 
must be fixed as In equation (4). Values of Nuq and Bsq are thus 
obtained for given values of x/d and p. 


Unfortunately, P varies along the passage. It will later be seen 
that qj ' is more convenient for use as a heat -flux parameter. 


2857 



CZ-2 back 


MCA TN 3016 


11 


The temperature and velocity distributions from figures 1 and 2 can 
be used for these q^uantlties except for the region outside the thermal 
boundary layer. The temperature distribution is^ of course, uniform 
outside the thermal boundary layer (f*" = . When the velocity distri- 

bution in that region is calculated for the case in which the velocity 
distribution is fully developed, it should be remembered that the 
properties are constant and have the values at ^h"^- values of 

5^^ >26, the velocity distribution curves outside the thermal boundary 

on the semllogarlthmic plot are straight lines having the slope of the 
curves at 

Nusselt numbers and Eeynolds numbers with properties based on tem- 
peratures other than the wall temperature can be obtained by use of the 
relations 


HAo = 

^iAo = 

PbAo = i/d-ptb"") 

PbAo = ^bAo = 

and so forth, where the perfect gas law with constant static pressure 
across the tube, constant Prandtl number, and constant specific heat 
have been assumed as in reference 6. The exponent d is obtained from 
viscosity data and is found to be about 0.68 for air and most gases. 

The following relations for Eeynolds numbers then result: 


E&b = Rsq 


= Ee^ 


d+1 


(^bAo) 

('t^bAp)^ 

A^Ap) 


where the continuity relation P-j^ U|j = has been used. 


(18) 

(19) 

Similarly, 


NUb = 


P 


(^bAp)*' 


( 20 ) 
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and 




Nub 


Another useful parameter is S viiich can he calculated from 

p(tbAo) 

q ' = — 

^^b AiAo) 


( 21 ) 


(22) 


Equation (22) can te obtained from the definition of ' and the 
perfect gas law for constant pressure, “toAb = ^bAo^ ^ 

constant along the tube. The important property of is that 

it does not vary along the length of the tube. 

It is desired to calculate the variation of bulk temperature and wall 
temperature along the tube for uniform heat flux. It can be shown, by use 
of the expressions qQ = h(tQ-t-jj) and itDXqQ = Pg(jt/4)D^ u-j^ Cp(t^-tj^), 

that dimensionless bulk temperature and wall temperature parameters are 
given by 



Equations (23) and (24) can be used for flow between flat plates as 
well as for tubes. 

The average Kusselt number for a given x/D Is calculated from 


^^i,av - px/D 

/ d(X/D) 

Jo ^*^1 


( 25 ) 


This equation is consistent vith the definition 
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r 


(to”t^)dX 


( 26 ) 


which la the usual way of 
and hulk temperatures for 
eq.uatlon (26) Is replaced 


(^0 - ■th)av 


defining the average difference between wall 
uniform heat flux. For uniform wall temperature 

[(^0-th)/(*loAoavj]^^ 

X 



which gives a Nusaelt numher consistent with equation (25) . For points 
heyond the entrance region the Nusaelt numbers and temperature distri- 
butions are calculated as outlined in appendix C. 


Development of Flow Boundary Layer 


Momentum equations for flow boundary layers . - The momentum equa- 
tion for the fluid in a given space can be written for steady flow in 
vector notation as 


F = 



v(pv • dA) 


(27) 


or the force acting on the fluid in a space is equal to the rate of 
flow of momentum out of the space (ref. 10, p. 233). When the X com- 
ponent of equation (27) is written for the previously illustrated 
differential annulus (6j^ replaced hy 5), there results 


dFx 



pu^ 


2irr dy) 



where the differential of the second Integral is the radial mass flow 
into the annulus at 5. The force dFx is composed of pressure forces 

on planes 1 and 2 and a shear force at the wall. There is no shear 
force at 6 because the velocity gradient is, by definition of the 
bound€iry layer, zero at the edge of the boundary layer. The force dF^ 

is then equal to 

- - Jt(rQ-S)^dp - 2rt roT^o*^ 

Substituting for dF^ gives, for equation (28), 
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In the region outside the boundary layer, the flow is frictionless, so 
that 

dp = - Pu^ dug (50 ) 


Substituting equation (30) in equation (29), rearranging, and integrating 
between the tube entrance and X give 



where the boundary- layer thickness has been taken as zero at X = 0. 
Equation (31) is the desired equation for flow in tubes relating the 
boundary -layer thickness to distance along the tube. 


Equations corresponding to equations (29) and (3l) can be derived 
for flow between flat plates. These equations are 



( 33 ) 
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Equation (32) vas first obtained bj yon Karman and has been used exten- 
sively in the incompressible form for calculating the grovth of the 
boundary layer for fluids flowing over surfaces (ref. 7, p. 66 ). 


The value of 5 must be known in terms of U 5 to solve equa- 
tions (31) and (33). This relation will be obtained in the section 
Conservation of mass^ constant properties. 


Dimensionless form of momentum equations for constant properties . - 
Equation (31) for tubes can be ^wi^tten Iji ^imensio form for constant 

fluid properties as 


X 

D 


1 

4 


\J 




lie 

2 


(5+/ro+)(2-6+/ro+) 


U 5 + d(u 6 '‘' ro'*’) - 



(34) 


Equation (34) can be verified by substituting the definitions of the 
dimensionless quantities into the equation. Equation (34) can be put 
into a form more convenient for computation by using the definition for 
the differential of a product. Thus 


^0 


+ 






(35) 


Substitution of equation (35) into the last term of equation (34) 
(the substitution can be more easily made if eq. (34) is first differen- 
tiated with respect to X) and simplification give 
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Equation (36) relates X/D to for various values of Ee 

for tubes. The relation between u"^ and y"*" is obtained from the 

curve for 3 = 0 in figure 2. The value of ro"^ for a given Ee and 

6+ will be obtained in the next section. The corresponding equation 
for flow between flat plates is found from equation (33) to be 



d r, 


(ug"^ - u'^)u‘^ dy 


Conservation of mass, constant properties . - In order to solve 
equations (36) and (37), u^"^ must be known In terms of Ee. This 

relation can be obtained from the law of conservation of mass. For a 
constant-area duct and constant density, the law of conservation of mass 
gives 

ui = Ut (38) 

or the buUc velocity is constant and equal to the uniform initial 
velocity. For a tube, equation (38) can be written 

2 P^O 

u. = — u(r^ - y)dy 


u(r„ - y)dy + — ^ (i 


(39) 
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Equation (39) can te written In dimensionless form as 


2ug+(r + - 5"^)' 

+ /„ + _+>j_+ . ^ ' 0 


Ee = — 7 I u (r^ - y )dy + 


•0 


(40) 


10 


■0 


The equation corresponding to equation (40) for flow "between flat plates 
is 

' 6 + 

Ee = 4 I u'*' dy"*" + 4ug'''(rQ''' - &■*■) (41) 


Equation (40) or (41) gives the relation between Ee and S'*' or Ug"*" 
for a given value of Tq''". 


ITusselt numbers, Eeynolds numbers, and friction factors; constant 
fluid properties . - Nusselt numbers for a developing velocity distribu- 
tion as well as a developing temperature distribution can be calculated 
from equations (10) and (12) or (14) and (16) as before. If the assump- 
tion is made that the thermal- and flow-boundary- layer thicknesses are 
equal (^h''' “ S"*" for Pr = l)^, equation (12) for tubes can be simplified, 
and for constant fluid properties (p = O) becomes 



A similar equation can be obtained for flow between flat plates . 

Eeynolds numbers - for flow in tubes and between plates can be calcu- 
lated from equations (40) and (41) . For constant fluid properties the 
Eeynolds number is constant along the length of the tube. 

The pressure drop between the entrance and a given section, which is 
required to calculate the friction factor, is obtained by integrating 
equation (30) . 

2 2 
P Ug P u^ 

Pi - p = 2 2 

Or, in dimensionless form, 

p' = 2(u§+ ^ 


(43) 

(44) 


^This assumption is used only for determining ti^''’ for the case where 
the Initial velocity and temperature distributions are uniform. 



18 


NACA OT 3016 


for tutes, or 


P' = 8(^5^ ■ t 




(45) 


for flat plates. The pressure-drop parameter p' In these equations is 
defined as 

. - (Pi - (4S) 


P = 


The average friction factor hased on the pressure drop for a given X/D 
is defined hy 


av 


- ^(Pj ' 

2XPU^^ 


(47) 


or 


■av - 2 


(X/D)Ee 


(48) 


This equation includes momentum changes caused hy velocity profile 
development . 

The local friction factor hased on the local pressure gradient is 
defined hy 

D(dp/dXl (49) 


f = - 


2pu> 


or 


f = 


_ 1 dp'/d(X/D) 


Ee'= 


(50) 


When f is to he found from equation (50) , the curve for p’ against X/d 
can be differentiated; or the curves, calculated from equation (56) or (37), 
for Ug"*^ ^o"*^ against 'L/T) can he differentiated and dp '/d (X/D) cal- 
culated from equation (44). The latter procedure proved to he the more 
accurate . 

The local friction factor hased on the local shear stress at the 
wall and excluding momentum changes due to a developing velocity profile 
is defined hy 

2x, 




"O 


P Ut^ 


( 51 ) 
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or 



(52) 


and, since Ee is constant along the tube and equal to 




8ro+ 


(53) 


where the relation between Tq^ and Ee for a given 5 is obtained 
from equation (40) or (4i). 


EESULTS AND DISCUSSION 

Heat Transfer with Uniform Heat Flux, Uniform Initial 
Temperature Distribution, Fully Developed Velocity 
Distribution, and Variable Fluid Properties 
for Gases; Prandtl Number of 0.73 

Development of thermal boundary layer . - Figure 3 illustrates the 
development of the temperature profile in a tube for uniform heat flux 
and a fully developed velocity profile. The distance along the tube 
for various thermal -boundary- layer thicknesses is calculated from equa- 
tion (4) and the temperature distribution in figure 1. The temperature 
distribution in the thermal boundary layer is obtained from figure 1; 
the temperature outside the boundary layer is uniform, inasmuch as no 
heat penetrates that region. The Eeynolds numbers are calculated from 
equations (ll) and (13). Figure 3 Indicates that for a fixed tube 
diameter at a given distance from the entrance, the boundary- layer 
thickness is less for high Eeynolds numbers than for low ones, or the 
distance required for a given degree of temperature profile development 
is greater for high Eeynolds numbers than for low ones. 

Although the temperature distributions in figure 3 are for p = 0, 
or constant fluid properties, similar distributions are obtained for 
variable fluid properties with large amounts of heat transfer, as may 
be Inferred from figure 4. The variation of the thermal entrance 
length, defined as the length required for the thermal boundary layers 
from the opposite walla to meet at the center of the passage, with 
Eeynolds number and heat flux ^ is shown in figure 4. The quantity 

q^ ' was calculated from equation (22), The quantity q ' rather than 

P is used as a heat -flux parameter Inasmuch as it is constant along 
the length of the tube. For the range of q^^ ’ shown, tQ/t-^ varies 

from about 0.4 to 3. Figure 4 indicates that heat flux has but a 
alight effect on the thermal entrance length. 
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Local Nusselt numbers for tubes . - Local Nusselt numbers hD/k^ 
with properties evaluated at the inlet temperature are plotted against 
X/D for various values of Ee^ and q,j^ ' (heat flux) in figure 5. The 

heat-transfer coefficient in the husselt number is based on the differ- 
ence between wall and bulk temperature ( see appendix A) . For 
obtaining the curves Nusselt numbers were calculated from eq^uations (lO), 
(20), and (21); Reynolds numbers, from equations (ll), (18), and (iS); 
q^', from equation (22); and X/d, from equation (4) for various values 

of the parameters Interpolation was necessary for 

obtaining the variation of Nu^^ with X/d for constant values of Re^ 
and qjj^ ' . 


The curves in figure 5 indicate, as ml^t be expected, that the 
Nusselt numbers (heat -transfer coefficients) close to the entrance are 
very high in comparison with the fully developed values farther down the 
tube at the ends of the curves. These values are high because 
the theimal boundary layer is thin and the temperature gradients conse- 
quently are severe near the entrance. At x/D = 0 the boundary- layer 
thickness is zero, so that the heat-transfer coefficient is Infinite 

at the entrance. For uniform heat flux this means that the temperature 
difference tQ - t-j^ must be zero at X/D = 0; for a finite temperature 

difference at x/D = 0 the heat transfer per unit area could not be 
uniform, since it would be infinite at the entrance. 

Inspection of the curves in figure 5 shows that the Nusselt numbers 
very nearly reach their fully developed values long before the thermal 
boundary layers have met at the ends of the curves. For instance, for 
q^’ = 0 (constant fluid properties) and Re^ = 100,000, the Nusselt number 

at x/d of 9 is within 2 percent of its fully developed value at X/D 
of 18.8. This indicates that, for practical purposes, the thermal entrance 
lengths can be considered to be less than half the values given in 
figure 4 . 

Comparison of the curves for q^ ' =0, 0.004, and -0.0025 indicates 
a decrease in NUj^ at a given x/d for heat addition to the gas and 

an increase for heat extraction from the gas. For both heat addition 
and extraction, Nuj^ increases with X/D for large values of X/D. 

These trends are clearly shown in figure 6, where Nuj^ is plotted for 
qj^ ' = 0.004 and -0.0025 for large values of x/D. For q^' = 0 (constant 

properties) the Nusselt numbers are constant beyond the entrance length. 

The Nusselt numbers for values of X/D beyond the entrance region were 
calculated by the method described in appendix C. For heat extraction 
from the gas the curves are cut at the points shown because the wall 
temperature reaches absolute zero at those points. The portions of the 
curves for which the temperature is below the liquefaction temperature 
of the gas should, of course, be disregarded because the properties no 
longer have the assumed variations with temperature. The Increases in 
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Kuaselt numl^er vith X/D are caused by the variation of properties along 
the tube and across the tube* The increases are due principally to the 
variation of which decreases along the tube for both heat addition 

to and extraction from the gas. 

It is clear from figure 6 that fully developed heat transfer, in 
the sense that the heat-transfer coefficient becomes independent of X/d, 
cannot be obtained when the properties are variable. Fully developed heat 
transfer can, however, be obtained in the sense that the relation between 
the Nusselt number and Reynolds number becomes independent of X/D for 
sufficiently large values of X/D when the properties are evaluated at 
the proper reference temperature, as shown in figure 7. The curves in 
figure 7 were obtained by cross-plotting the curves for ’ = 0 from 

figure 4, They can be used for ^ 0 by evaluating the properties, 

including density, in the Nusselt number and Reynolds number at the 
reference temperatures t^ ^ " "^b^ "^b 1^ figure 8. The 

values of x were computed by selecting values of t^ such that the 
values of Nu^ and calculated from the curves in figure 5, from 

equations (18) to (2l), and from the definition of t^, fall on the lines 
for q^ ' = 0 in figure 7 . The curves in figure 8 indicate that the 

reference .temperature for points close to the entrance is slightly 
closer to the bulk temperature then the reference temperature for fully 
developed heat transfer, although its variation with x/D is not large. 

Effect of variation of heat transfer across boundary layer on 
Nusselt number for small 5^"'" and constant fluid properties . - For large 

values of {tq^ for fully developed flow), the effect of the radial 

variation of heat transfer on temperature distributions was checked in 
reference 2 and found to be small. However, for low values of 
such as occur very close to the entrance, the flow is nearly laminar, 
so that the temperature profile in the boundary layer is not flat as it 
is for higher values of The temperature distribution in the 

thermal boundary layer for laminar heat transfer from a flat plate can 
be closely represented by a cubic parabola (ref. 7, p. 89), which means 
that the heat transfer per unit area varies as 1 - y^/6^^, as shown in 

appendix B. This variation of heat transfer is used herein for turbulent 
heat transfer for small 6^“*“ in order to obtain the effect of variation 

of on the Nusselt numbers in the entrance region. The expression 

for the temperature distribution with this variation of q/qQ is obtained 
in appendix B. 
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The variation of Nu with x/D for variable heat transfer per 
unit area across the boundary layer (i/Iq = 1 - shown in fig- 

ure 9 For comparison, the curves from figure 5(a) for constant heat 
transfer (q./q.Q = l) are shown by dotted lines. The curves in figure 9 

indicate that the variation of heat transfer per unit area across the 
boundary layer has a negligible effect on the relation between Kusselt 
number and X/d, so that the use of the temperature distributions in 
figure 1 appears Justified, even for low values of • 

Wall and bulk temperat ure distributions for tubes . - Wall and^ 
bulk temperatures along a tube with uniform heat transfer as calculated 
from eq.uations (23) and (24) are plotted in dimensionless form in 
figure 10. As mentioned previously, the wall and bulk temperatures are 
eq,ual at the entrance for uniform heat transfer. The bulk temperature 
varies linearly with X/D Inasmuch as the specific heat is assumed 
constant . For ' = 0 the wall and bulk temperature curves are 

parallel beyond the entrance region. For both heat addition 

(q^’ = 0.004) and heat extraction (ij^' = -0.0025) from the gas, the 

wall £ind bulk temperature curves converge beyond the entrance region 
because of the increase in heat-transfer coefficient with X/D for 
those portions of the curves. 

Average Husselt numbers for tubes . - The variation of average 
Nusseit number Nu^ calculated from equation (25) is 

shown in figure 11.^ The trends are similar to those for local Nusseit 
numbers but the changes are more gradual. 

The relation between N'^x,av ^®x,av various values of 

x/d is given in figure 12. The curves were obtained by cross -plotting 
the curves for q^ ’ = 0 in figure 11(a) and can be used for qj ’ / 0 

by evaluating the properties in the Nusseit number and Eeynolds number 
at . i(to^a, - where I la given In figure 13. 

The values of x were computed by selecting values of tx,av that 

the values of Nu^ ^v ®®x,av calculated as outlined in appendix C 

fall on the lines for q^ ' =0 ’in figure 12. Figure 13 indicates that 

the value of x increases continuously with X/D for both heat addition 
to and extraction from the gas. These trends are Indicated experimentally 
for heat addition in reference 11 . where a constant reference temperature 
used for a range of values of X/D was found to give average Nusseit 
numbers for hi^ heat fluxes at low values of x/d which were slightly 
jjj^gher than those for low heat fluxes. The values of x given in 
figure 13 indicate that for large values of X/D, x>1, so that for heat 
addition t., „„ is higher than the average wall temperature although 

not necessarily higher than the exit wall temperature. The variation 
of X with X/d is caused by the variation of properties along the 
tube. A more satisfactory correlation would perhaps be obtained by 
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using X only to account for the variation of properties across the 
tuhe and using another factor to account for the variation along the 
tube. The curves in figures 12 'and 13 mi^t, however, provide a con- 
venient means of utilizing the results of the analysis for some purposes. 

Nusselt numhera and temperature distrlhutions for flow between flat 
plates . - Nusselt numbers, wall temperature distributions, and bulk tem- 
perature distributions for flow between flat plates with uniform heat 
transfer and fully developed velocity distribution are plotted In 
figures 14 to 16. The curves agree closely with the corresponding curves 
for flow in a tube. For instance, the curves for local Nusselt numbers 
plotted in figure 14 are between 3^ and 5 percent higher in the fully 

developed portions of the curves than the corresponding curves for tubes 
in figure 5. The portions of the curves close to the entrance agree 
even more closely. This agreement indicates that results for turbulent 
heat transfer in tubes can be used, with small error, for flat plates 
when an eq^uivalent diameter eq.ual to twice the plate spacing is used. 

It should be emphasized that in applying the results of the fore- 
going analysis the same assumptions regarding the variation of properties 
with temperature should be made as wep^e made in the analysis, that is, 
constant specific heat, and viscosity and thermal conductivity both 
proportional to t'^*^®. 

It should also be mentioned that the results in this analysis do 
not include the low Peclet number effect described in reference 8 inas- 
much as the temperature distributions in figure 1 do not Include that 
effect. The Nusselt numbers at He^ = 10,000 may therefore be slightly 

high, but the ratios of local to fully developed Nusselt numbers should 
be accurate . 

Effect of frictional heating on entrance length . - The effect of 
frictional heating on the X/D for a given can be approximately 

accounted for by replacing tg+ and t+ by T5+ and T+, respectively, 
in equation (4) and setting p/Pq = l/(l - PT+ - ou+^) (ref. 9). In 

using that procedure it is assumed that the total-temperature gradient is 
zero at the edge of the thermal boundary layer, a good assumption for 
Prandtl numbers close to unity. 

The variation of thermal entrance length (X/D)^q with Eeynolds 

number and heat flux with frictional heating (a = 0.00025) is shown 
in figure 17. Mach numbers for air at the tube center at (^/®)he 

Indicated in the figure. The Mach numbers were calculated from 
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which can he obtained from the definitions of M, T^,, and a. Com 

parison of figure 17 with figure 4 indicates that frictional heating 
increases the entrance length, although the Increase is not large for 
subsonic flow. Inasmuch as Nusselt numbers for fully developed heat 
transfer are not affected by frictional heating (ref. 9), the effect 
of frictional heating on the relations between Kusselt number and X/D 
probably is slight when the total rather than the static bulk tempera- 
ture is used in the definition of the heat-transfer coefficient. 


Heat Transfer with Uniform Wall Temperature, Uniform Initial 
Temperature Distribution, Fully Developed Velocity 
Distribution, and Constant Fluid Properties for 
Gases; Prandtl Number of 0.73 

Relations among Nusselt number, Reynolds number, and X/D for the 
case of tubes, as calculated from equations (9) to (11), are- plotted in 
figure 18. Comparison of figure 18 with figure 5 indicates that the 
Nusselt numbers for uniform wall temperature are very slightly lower 
than those for uniform heat flux. The heat transfer per unit area qg 

must be Infinite at X/D = 0 for uniform wall temperature. This occurs 
only for an infinitesimal distance dX, however, so that a finite amount 
of heat is transferred . 


Heat Transfer with Uniform Heat Flux, Uniform Initial Velocity and 
Temperature Distributions^ and Constant Fluid Properties 
for Gases; Prandtl Number of 0.73 

For calculating the variation of Nusselt number with x/d for 
t^is case the flow-boundary-layer thickness as well as the thermal- 
boundary-layer thickness should be known for various values of X/D. 
However, in the calculation of Nusselt numbers, it is considered 
sufficiently accurate to assume the thicknesses of the thermal and 
flow boundary layers to be equal. Inasmuch as the Prandtl number is 
close to unity. 

Relations among Nusselt number, Reynolds number, and x/D for 
tubes, as calculated from equations (lO), (42), (40), and (4), are 
plotted in figure 19. Comparison of figure 19 with figure 5 indicates 
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that the Nusselt numhers for a uniform Initial velocity distribution are 
higher than those for a fully developed velocity distribution. This 
result might be expected because of the higher friction in the case of 
the uniform initial velocity profile. 

Average Nusselt numbers as calculated from equation (25) are plotted 
in figure 20. Figure 20(b), which is a logarithmic plot, compares the 
predicted relations between Nu^^ and X/D with the empirical relation 

Nu a(x/D)*^-^. Curves having a slope of -0.1 on the log-log plot agree 
closely with the predicted lines for X/D between 6 and 60. The agree- 
ment of the data given in. reference 11 with the empirical relation, for 
values of X/d as high as 120 might be due to the use of a constant 
rather than a variable reference temperature as given in figure 13, 

Nusselt numbers for flow between flat plates are given in figure 21, 
Comparison of figure 21 with figure 19 indicates very close agreement 
between the Nusselt numbers for tubes and flat plates. 


Heat Transfer with Uniform Wall Temperature, Uniform Initial 
Velocity and Temperature Distributions, and Constant Fluid 
Properties for Gases; Prandtl Number of 0.73; 
and Comparison with Experimental Data 

The variation of Nusselt number with X/d for this case is given 
in figure 22. The curves agree closely with those for the corres- 
ponding case for uniform heat flux given in figure 19. The Nusselt 
numbers for uniform wall temperature are slightly lower than those for 
uniform heat flux. 

A comparison between analytical and experimental results for uni- 
form wall temperature and uniform initial velocity and temperature dis- 
tributions is given in figure 23. The experimental curves represent 
data from reference 12 for air flowing In a tube at uniform temperature 
with a bellmouth entrance and a screen at the entrance and should repre- 
sent the case analyzed in this section. The screen at the entrance 
should Insure a turbulent boundary layer throughout the tube. The 
figure shows that the agreement "between analytical and experimental 
results for this case is very good. 

The lack of a comparison of the present analytical results for a 
uniform wall temperature and a fully developed velocity distribution 
with the data from reference 12, which were taken with a long approach 
section at the entrance, should perhaps be explained. It appears that 
the sharp edge at the beginning of the long approach section produced 
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large disturbances which were not damped out in the length of approach 
section used. The data in reference 12 for heat-transfer coefficients 
with a sharp entrance without an approach section Indicate such results. 
For that case the heat-transfer coefficients are hl^er throughout the 
tube than those obtained with a bellmouth entrance. 

Average Nusselt numbers calculated from eq.uation (25) are plotted 
in figure 24. As in the case of local Fusselt numbers, the average 
Nusselt numbers for a uniform wall temperature in figure 24 are slightly 
lower than those for a uniform heat flux in figure 20. 


Heat Transfer with Uniform Heat Flux, Uniform Initial 
Temperature Distribution, Fully Developed Velocity 
Distribution, and Constant Fluid Properties for 
Llqtuld Metals; Prandtl Number of 0.01 

Temperature distributions for fully developed flow calculated as 
described in reference 8 and plotted in figure 25 were used in the 
thermal boundary layer for calculating the variation of Nusselt number 
with X/D for liquid metals. The velocity distributions were obtained 
from the curve for P = 0 in figure 2. The curves for Nu against 
x/D in figure 26 were calculated from equations (4) and (lO) to (12) 

(P = O). An unexpected feature of these curves is the sli^t Increase 
of Nusselt number with X/D at large values of X/D and Eeynolds num- 
ber. This increase is due to the change in shape of the temperature 
profile with for low Prandtl numbers. Thus, as bjj'*' increases, 

the temperature gradient at the wall increases slightly in some 
instances because of the added turbulence in the boundary layer at 
hi^er values of 6^^+. This effect is indicated experimentally by 

the data in reference 13. 

A plot of Nusselt number against Pe/(X/D) for various Peclet 
numbers is given in figure 27. If the heat transfer contributed by 
eddy diffusion were neglected, the curves for various Reynolds numbers 
would fall essentially on a single line as in reference 5. Examination 
of the curves in figure 27 indicates that for Reynolds numbers below 
1000 the curves could be represented approximately by a single line. 
This line is slightly higher than the one obtained in reference 5 inas- 
much as a uniform wall temperature rather than a uniform heat flux was 
postulated in that reference. 
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Friction Factors with Uniform Initial Velocity Distribution 
and Constant Fluid Properties 

Three friction factors -are discussed In this section: f^^, based 

on the shear stress at the wall; f, baaed on the static-pressure gradient; 
and fg^^ based on the pressure drop along the tube. The friction factor 

fx is calculated from equation (53); f, from equation (50); and fay^ 

from equation (48). The Refolds numbers and values of x/d corres- 
ponding to these friction factors are calculated from equations (41) and 
(36) , The velocity distributions in the flow boiundary layer are 
obtained from the curve for 3=0 in figure 2 . 

The friction factors for a tube based on the shear stress at the 
wall and on the static-pressure gradient are plotted in figures 28 and 
2S, respectively. The friction factors based on the pressure gradient 
are higher than those based on the shear stress in the entrance region 
because the former Include the pressure decrease caused by the increase 
in momentum of the fluid as the velocity profile develops as well as 
the pressure decrease caused by the shear stress at the wall. For 
fully developed flow the two friction factors are identical. Average 
friction factors based on the pressure drop are plotted in figure 30 
and are, of course, higher than the local values. 

The computed values of f and are subject to some Inaccuracy 

because of errors associated with measuring the slope of a curve (eq. 

(50)) and because the two terms in equation (44), which are subtracted, 
are nearly equal close to the entrance. The results presented should, 
however, be well within the limits of accuracy of experimentally deter- 
mined friction factors. 

Friction factors for flow between flat plates are plotted in 
figures 31 to 33 and are found to be very similar to those for a tube. 


Comparison of Local Uusselt Numbers and Friction Factors 

A comparison of local Nusselt numbers and friction factors for 
various cases 1s presented in figure 34, where the ratio of local 
Nusselt number or friction factor to the fully developed value is 
plotted against X/D for a Reynolds number of 100,000. The quantity 
Nu/Nu^ is smaller than f/f^ for a given x/d in all cases except 

that of heat transfer to a liquid metal (Pr = O.Ol). The larger values 
of Nu/Nu^ for a liquid metal are apparently caused by the fact that 

the temperature profile at low Prandtl numbers is similar to that for 
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laminar flow. At very high Reynolds numbers, where the dip occurs in 
the curves for Wusselt number plotted against X/D, and at very low 
Reynolds numbers, the values of Nu/Nu^j may be lower than f/f^ for 

liquid metals (fig. 26). For gases the values of Ku/Nu^^ are lower 

than f/% for a given value of X/D because of the momentum pressure 
loss for velocity profile development included in the friction factor 
f; there is no counterpart for this momentum pressure loss in heat 
transfer . 

It should be mentioned that the preceding conclusions may apply 
only to the case analyzed, that is, to the case where the boundary layer 
is turbulent throughout the tube . The turbulent boundary layer might be 
produced by a bellmouth with a strip of sand paper at the entrance or 
by screens. For flow in which the boundary layer is partially laminar, 
or for which there are large disturbances at the entrance, such as might 
be caused by a right -angle-edge entrance, the conclusions mi^t be 
altered . 


SUMMARY OF RESULTS 

The following results were obtained from the analytical investiga- 
tion of the entrance region for turbulent heat transfer and flow in 
smooth tubes and between parallel flat plates: 

1. The thin thermal boundary layers and consequently severe tem- 
perature gradients at the wall near the entrance produced hl^ heat 
transfer near the entrance. Approximately fully developed heat transfer 
was, in general, attained in an entrance length less than 10 diameters. 

2. The effect of variable properties on the local Nusselt number 
correlation for gases with fully developed velocity profiles was to 
decrease the Russelt number with conductivity based on the inlet tem- 
perature for heat addition to the gas and to increase it for heat 
extraction from the gas at a given ratio of distance from the entrance 
to tube diameter X/D and inlet Reynolds number. For large values of 
X/D, the Nusselt number with conductivity based on the inlet temperature 
increased continuously with x/D for both heat addition to and extrac- 
tion from the gas. A heat-transfer coefficient which is independent of 
X/D Is therefore not obtained for heat transfer with variable properties, 
even in the fully developed region. 

3. The reference temperature used for evaluating the properties 
in order to eliminate the effect of heat flux on the curves for local 
Nusselt number against Reynolds number at various values of X/D varied 
only slightly with x/d. The reference temperature for average Nusselt 
numbers, however. Increased continuously with X/D for heat addition 
to the gas and was higher than the average wall temperature but lower 
than the exit wall temperature at very large values of x/D. 
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4. Th© variation of Hfusselt number with X/D for a gas with uni- 
form initial temperature and velocity distributions and uniform heat flux 
agreed with the experimentally determined relation, Nu«c for 
values of X/D between 6 and 60, 

5. The predicted variation of Nusselt number with X/d for air 
with a uniform wall temperature and uniform initial velocity and tem- 
perature distributions agreed very well with the experimental values 
from reference 12 , 

6. The variation of Nusselt number with Graetz number for liq^uid 
metals (Prandtl number of 0,01) was represented approximately by a single 
line for Pec let numbers below 1000. 

7. The ratio of friction factor based on the static- pres sure 
gradient in the entrance region to the fully developed value was greater 
than the corresponding ratio of Nusselt numbers for gases in the entrance 
region because of the momentum pressure loss associated with the velocity 
profile development, 

8. Essentially the same variatlohs of Nusselt number and friction 
factor with X/D were obtained for flow between parallel flat plates 
as were obtained for a tube when the equivalent diameter equal to twice 
th© plate spacing was used for th© plates. 

Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, July 21, 1953 
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APPENDIX A 
SYMBOI^ 


The following syrnhole are used 

A area, sq ft 

Cp specific heat of fluid at 

c specific heat of fluid at 

^0 Btu/(lh)(°F) 


in this report; 

c onstant pre s sure , Btu/ ( Ih ) ( °F ) 
constant pressure at wall. 


D Inside diameter of tube, ft 

d exponent; value depends on variation of viscosity of fluid 

with temperature 

g acceleration due to gravity, 32.2 ft/sec^ 

h local heat-transfer coefficient, qQ/tQ - t^, Btu/(sec)(sq ft)(°F) 


h^^ average heat-transfer coefficient, Q.o^av/^^0 ~ ^^av 

J mechanical equivalent of heat, 778 ft-lh/Btu 

h thermal conductivity of fluid, Btu/(sec)(sq ft)(°F/ft) 

thermal conductivity of fluid evaluated at tr^^, 
Btu/(sec)(sq ft)(°F/ft) 

kj^ thermal conductivity of fluid evaluated at tj^, 

Btu/(sec)(sq ft)(°F/ft) 

kQ thermal conductivity of fluid evaluated at tQ, 

Btu/{sec)(sq ft)(°F/ft) 

kj thermal conductivity of fluid evaluated at t^, 

Btu/(sec)(sq ft)(°F/ft) 

k^ thermal conductivity of fluid evaluated at 

Btu/(sec)(sq ft)(°F/ft) 
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M Mach numher at tube center at X/D 

n constant 

p static pressure 

static pressure at inlet^ Ih/sq. ft ahs 


^0 


^0, 


av 


^0 

T 


^h,av 


rate of heat transfer toward tube center per unit area, 

Btu/ (sec) (sq. ft) 

rate of heat transfer at inside wall toward tube center per 
unit area, Btu/(sec) (sq. ft) 

average rate of heat transfer at inside wall toward tube 
center per unit area for a given X/D, Btu/(sec)(sq ft) 

radius, distance from tube center, ft 

inside tube radius or one-half distance between plates, ft 

absolute total temperature, °R 

absolute total temperature at tube center, 

total temperature outside thermal boundary layer, 
absolute static temperature, ^ 

bulk static temperature of fluid at cross section of tube, ^ 

average bulk temperature, °R 

inlet static temperature of fluid, ^R 


t^ reference temperature for local Nusselt and Reynolds numbers, 

x(tQ - t^) + °E 

t„ reference temperature for average Nusselt and Reynolds numbers, 

^ £1 V . - 

^^^0,av “ ^,av' ^b,av 

t^ temperature of fluid outside thermal boundary layer 
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tQ absolute wall temperature, 

■to,av average wall temperature, °E 

(to-t"b)gy average difference tetween wall and tulk temperature, 
u time-average velocity parallel to axis of channel, ft/sec 

hulk velocity at cross section of tube, ft/aec 

Uj velocity of fluid at inlet, ft/sec 

u_ velocity outside flow boundary layer, ft/sec 


X 

X 

7 

6 

^h 

e 

h 

^^b 

^^i 

^'x 

^^x,av 

P 


velocity vector, ft/sec 
distance from entrance, ft 

number used in evaluating arbitrary temperature in tube, 
^x ^x,av 

distance from wall, ft 

flow-boundary -layer thickness, ft 

thermal-boundary- layer thickness, ft 


coefficient of eddy diffuslvity, sq ft/sec 
absolute viscosity of fluid, (lb)(sec)/sq ft 

absolute viscosity of fluid evaluated at (lb)(sec)/sq ft 

absolute viscosity of fluid evaluated at t^, (lb)(sec)/sq ft 

absolute viscosity of fluid evaluated at t^^-, (lb)(sec)/sq ft 

absolute viscosity of fluid evaluated at 

(lb)(sec)/sq ft 

absolute viscosity of fluid evaluated at tQ, (lb)(sec)/sq ft 
mass density of fluid, (lb) (sec^)/ft^ 
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mass density of fluid evaluated at t^j, (lb) (sec^ )/ft^ 
Pb,av mass density of fluid evaluated at t^ (lb) (sec2)/ft'^ 

Pj^ mass density of fluid evaluated at (lb) (sec^)/ft^ 

P^^ mass density of fluid evaluated at t^.^ (lb) (sec^)/ft^ 

P^ av mass density of fluid evaluated at t^j. (lb) (sec^)/ft^ 

Pq mass density of fluid evaluated at tQ, (lb) (sec^)/ft^ 

T shear stress in fluid, Ib/sq ft 

Tq shear stress in fluid at wall, Ib/sq ft 

Dimensionless groups: 

f friction factor based on static-pressure gradient, 

2p u^,'^ 
D(Pi - p) 

fg^ friction factor based on static-pressure drop, 

2 X p u^'^ 

f^ fully developed friction factor 

friction factor tased on shear stress at wall, 2 Tq/(p 

Gz Graetz numher, Re Pr/ (X/D) 

Nu Russelt nuznher, h D/k 

average Nusselt numher for constant properties defined hy 
eq. (25) 

Nu-j^ Nusselt numher with thermal conductivity evaluated at t-^, 

h D/k^ 

Nu^ fully developed Nusselt number 

NUj_ Nusselt numher with conductivity evaluated at t^ 
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Nu 


i,av 


Nu, 


Nu, 


Nur 


P’ 


avarage Nusssl't nunibsr defined by eq. (25) with conduc- 
tivity evaluated at tj^ 

Nusselt number with conductivity evaluated at t^j. 
average Nusselt number with conductivity evaluated at 
Nusselt number with thermal conductivity evaluated at Iq 
pressure-drop parameter, (pj^ - p)pD'^/n^ 


Pe 

Pr 

1i' 

Re 

Re-jj 

Rsj^ 

Rej 

^®x,av 

Reo 


Peclet numlDer; Ee Pr 
Prandtl number, Cp g ti/k 

Prandtl number with properties evaluated at tg 
heat-transfer parameter, <lo/(®p 6 P^u^t^) 

Reynolds number, P u-|j I>/\i 

Reynolds number with density and viscosity evaluated at t-jj, 
Reynolds number with density and viscosity evaluated at t^^, 

PiV 

Pi 

Reynolds number with density and viscosity evaluated at t^j., 

Px% Vpx 

Reynolds number with density and viscosity evaluated at t^^^av^ 

E ^ Jd_ 

Pb,av Px,av 

Reynolds number with density and viscosity evaluated at tg, 

Po% D/Po , 

* 

tube radius parameter, 7 - — rg 
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ro 

CD 

ui 




Uv 




U5-^ 


total-temperature parameter, 


(tp - T)Cp g Tq ^ 1 - T/tp 

<1oVVpo ^ 


i 1 - - 

P V t 


temperature parameter, 


(tp - t)Cp g Tp ^ 1 - t/tp 

^oV'^'o/^o ^ 


tulk-temperature parameter 

i(‘-a 


•K‘-3 


velocity parameter, 


u 


V'^o/^o 


bulk-velocity parameter. 


% 


V'o/Po 


value of u'*' at tube center 




0/^0 


he 


a 

P 

6+ 


thermal entrance length divided by diameter 


-^"^/^/Pp 

wall distance parameter, ^ — - y 

^^o/Po 

compressibility or frictional heating parameter, '^p/(2eJCptpPp) 
heat-transfer parameter, <lp V”^o/^0 /^^p ^ "^0 ^0^ 

a/V^ 

dimensionless flow-boundary- layer thickness, ^ 


dimensionless thermal -boundary-layer thickness. 


toT^ 


An Integral sign having open brackets for the upper limit 

indicates that the variable of integration should be used as the 
upper limit . 
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APPENDIX B 


ANALYSIS INCLUDING EFFECT OF VAPIATION OF HEAT TE^SFER 
ACROSS BOUNDARY LAYER FOE SMALL WITH 

CONSTANT FLUID PROPERTIES 

As was msntioned in 'tlia section Effect of variation of heat transfer 
across Boundary layer on Nusselt numBer for small Bji'*' and constant 
fluid properties^ the effect of variation of t/q.q across the Boundary 
layer for small Bji"*” should Be checked . For laminar flow over a flat 
plate the temperature distriBution in the thermal Boundary layer can Be 
closely approximated By a cuBic paraBola. For a cuBic paraBola^ 

= 1 - as shown in the following: 


For laminar flow, 

q, = - k dt/dy = qo(^ " 


Integration of equation (Cl) gives 


- (^^0 - t ) 


36 


^0 


(Cl) 


(C2) 


Evaluation of equation (C2) at the edge of the Boundary layer gives 


^ (to - ts) = I 6j^ 

ComBlnation of equations (C2) and (C3) results in 

_ 3 . i /xf 

to - tg 2 6jj 2 \^Y^j 


(C3) 


(C4) 


which is the same as the equation on page 89 in reference 7. The 
assumption that <l/<lo = 1 ■ /^h^ therefore leads to a cuBic paraBola 

for the temperature distriBution for laminar flow. The same assumption 
is used herein for turBulent flow for small The equation for 

turBulent heat transfer can be written 
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q. = - (k + Cp g pe) dt/dy (C5) 

With the assumption for q/q.Q and e = n^uy for flow close to the wall 
(ref. 6), equation (C5) becomes, in dimensionless foim. 


1 






(C6) 


, r 

where t"^ represents the temperature distribution with 
2 2 

a/io = 1 - /^h"^ as distinguished from which is the dimensionless 

temperature parameter for i/<Lq = 1. Separation of variables and 
integration of equation (C6) result in 


t 


+ ' 



llL 


Pr 


+ n^u''‘y"'' 




+2 ^ + 
y 

^ + n^u+y'*’ 


(C7) 


The first term on the right side of equation (C7) is t'*’ for constant 
properties (ref. 6), so that 


t 


+ ' 




y+^ dy~*~ 

-i- + n^u"*’y* 
Pr 


(C8) 


The values of t"*" are obtained from figure 1. For obtaining the 
Nusselt numbers the values of from equation (C8) are used in 

equations (12) and (4) in place of f*". 
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APPENDIX C 


CALCULATION OF HUSSELT NUMBERS, REYNOLDS NUMBERS, AND TEMPERATURE 

DISTRIBUTIONS FOR FLOW BEYOND ENTRANCE REGION 

For points teyond the entrance region the relation between Nusselt 
number and Reynolds number with the properties, Including density, 
evaluated at the reference temperature t^^ = i(to - t^j) + t^, where 
X = 0.4 for heat addition and 0.6 for heat extraction from the gas, 
can be obtained from the curve for fully developed flow in figure 7. 
From equations (18) to (21) and the definition of t^ it can be shown 

that 



where d = 0.68 for air and most other gases. The value of Nuj^ for 
a given value of Re^ can be obtained from equations (Bl) and (B2) 
provided the values of t^^/t^ and to/t^ are known. The quantities 
t-^j/t^ and tQ/t^ can be written in tenns of the quantities plotted in 
figure 10 as 



(see eq. (23)) and 

■tb 

In the procedure for obtaining the value of ^qJ\> quantity in 

brackets in the numerator is first determined by extrapolating a known 
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J 


curye for that quantity plotted against X/D. The value of Nuj_ can 
then he obtained from equation (B2). With that value of Nuj_ a new and 


more 

tion 


accurate value of 

(24) . The calculation of NUj^ 


— ^ can he calculated from equa- 
can then he repeated using the nev 


value of 


[(to-ti)/ti] ^ 


Average Nusselt numbers with the conductivity evaluated at the 
inlet temperature can be calculated from equation (25), where 

the values of Nuj^ are calculated as outlined previously. It is 

desirable to calculate average Wusselt numbers and Heynolds numbers with 
properties evaluated at the average bulk or the average wall temperature 
for given values of x/D. The average bulk temperature rise for a given 
x/d is one-half the temperature rise at X/d, since the specific heat 
is assumed constant. Therefore, 



From the definitions of hg^^, NUj^ Ee^, and Pr, it can be shoxm that 


^0,av _ ^h.av ^ 


(B6) 


Average Nusselt numhera and Reynolds numhers with properties evaluated 
at t-jj tQ or any arbitrary temperature 

^x,av ~ ^Ao^av “ ^h,av^ ^h,av 
can he calculated hy use of equations (B5) and (B6) . 
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CZ-6 back 



FlfOire 3. - Development thermal boundary layer In entrance region of tube for a gae. Pr, 0.73j 0. Uniform heat flux 

uniform initial temperature distribution, and fully developed velocity distribution. 





Figure 4. - Variation of thermal entrance length with Reynolds number and heat flux for gas flowing In a tube, 
form heat flux, uniform Initial temperature distribution, and fully developed velocity distribution. 






2857 


NACA TN 3016 


45 



Figure 5. - Variation of Nusaelt number Nu^ vlth x/d and Reynolds number for gas flowing in a tube. Pr, 0,73. Uniform heat flux, 
uniform initial temperature distribution, and fully developed velocity distribution. 
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Figure 5. - Concluded. Variation of Nusseit number Nuj^ with X/D and Reynolds number for gas flowing in 

a tube. Pr, 0.73. Uniform heat flux, uniform initial temperature distribution, and fully developed 
velocity distritiution . 
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X 


Figure 8, - Variation of x, defined "by = x(tQ - t^) + with X/D for use with figure 
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Figure 10. - Concluded. Variation of dimensionless wall and bulk temperatures with x/d and 
Reynolds number for gas flowing in a tube. Pr, 0.73. Uniform heat flux, uniform initial 
temperature distribution, and fully developed velocity distribution. 
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temperature distrltutlon, and fully developed velocity distribution. 
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x/b 

(a) » 0 (constant properties). 


Figure 14. - Variation of Nusselt number with X/D and Reynolds number for gas flowing between parallel 
flat plates. Pr, 0.73. Uniform heat flux, uniform Initial temperature distribution, arid fully devel- 
oped velocity distribution. 
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(b) » 0.004 (heat addition to gas). 

Figure 14. - Continued. Variation of Nusselt number with X/D and Reynolds number for gas flowing 
between parallel flat plates. Pr, 0.73. Uniform heat flux, uniform Initial temperature distribution, 
and fully developed velocity distribution. 
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x/d 

(c) = -0.0025 (heat extraction from gas). 


Figure 14. - Concluded. Variation of Nusselt number with X/D and Reynolds number 
for gas flowing "between parallel flat plates. Pr, 0.73. Uniform heat flux, uni- 
form Initial temperature distribution, and fully developed velocity distribution. 
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Figure 15, - Variation of dlroenalonlega wall and bulk temperatures with X/D and Re, 7 nolda number for gas 
flowing between parallel flat plates. Pr, 0.73. IMlforra Initial temperature distribution and fully 
developed velocity diatrlbution . 
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x/d 

(a) q^' » 0 (constant properties). 

Figure 16. - Variation of average Nusselt numljer defined by equation (25) with X/D and Reynolds number 
for gas flowing between parallel flat plates. Pr, 0-73. Uniform heat flux, uniform initial tempera- 
ture distribution, and fully developed velocity distribution. 
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Reynolds number for gas flowing between parallel flat plates. Pr, 0.73. Uniform heat flux, uni- 
form initial temperature distribution, and fully developed velocity distributions. 
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(c) = -0.0025 (heat extraction from gas). 

Figure 16. - Concluded. Variation of average Nusselt number defined by equation 
(25) with x/D and Reynolds number for gas flowing between parallel flat plates. 
Pr, 0.73. Uniform heat flux, uniform initial temperature distribution, and 
fully developed velocity distribution. 
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Figure 18. - VfiLrl&'tlon of Ifussel't nuoiber vltih x/D and Reynolds number for gas flowing In a 
tube. Pr, 0.73. Uniform wall temperature, uniform initial temperature distribution, lully 
develox>ed velocity distribution, and constant properties. 
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gure 19. - Variation of Husselt number vlth x/d and Reynolds number for gas flowing in a tube 
uniform Initial temperature and velocity distributions, and constant properties. 
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(a) Rectangular Cartesian coordinates. 

Figure 20. - Variation of average Nusselt number CEilculated from equation (25) with X/d and Reynolds nvuflaber for gas flowing 
in a tube. Pr, 0.73. Uniform heat flux, uniform Initial temperature and velocity distributions, and constant properties. 
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gas ^Inwing in a tuhe. Pr j 0.73. Unl^^or:^ heat ^lux, uriforin Initial temperature and velocity distributions, and constant 
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Figure 21. - Variation of T^sselt number with X/D and Reynolds number for gas floving 
between parallel flat plates. Pr, 0.73. Uniform heat flux, uniform initial tempera- 
ture and velocity distributions, and constant propeirtles , 
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Figure 22. - Variation of Nusselt number with x/d and Reynolds number for gas flowing In a tube. Pr, 0.73. Uniform 
wall temperature, uniform Initial velocity and temperature distributions, and constant properties. 
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Figure 23. - Comparison of predicted and experimental results for air floving in a tube, Pr, 0.73. Uniform wall 
temperature, uniform initial velocity and temperature disti ibutions, and constant properties. 
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Figure 25. - Predicted generalized temperature distribution for flow of liquid metal with heat 
transfer. Data from re^'erence 8. Pr, 0.01. 
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Figure 26- - Variation of Nusselt number with X/D and Reynolds number for liquid metal flowing in 
tube. Prandtl number, 0.01. Uniform heat flux, uniform initial temperature distribution, fully 
developed velocity distribution, and constant properties. 



Figure 27. - Variation of Nusselt number with Graetz number and Peclet number for liquid metal flowing in a tube. Pr 0.01. 
Uniform heat flux, uniform initial temperature distribution, fully developed velocity distribution, and constant properties 
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Figure 29, - Variation of friction factor based on statlc-presaure gradient with X/D and Reynolds number for flow In a tube. 
Constant properties. 
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Figure 30. - Variation of average friction factor based on static-pressure drop with X/D for flow In a tube. Constant properties. 
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Figure 33. - Variation of average frict^ 
with x/d and Reynolds number for flc 
properties . 



1 factor based on static -pres sure drop 
between parallel flat plates. Constant 
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